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Cosmological Phase Transitions 

Jan Jalmar (Janne) Ignatius 
University of Helsinki, 1993 



Abstract 

It is generally believed that several phase transitions have taken place in the early 
Universe. The effects of cosmological phase transitions may well have been crucial for 
the evolution of the Universe, and thus for the existence of life as we know it. 

The cosmological phase transitions investigated here are related to strong and 
electroweak interactions. When the Universe was about 10~^^ seconds old and the 
horizon radius equaled one centimeter, symmetry between weak and electromagnetic 
interactions was broken. It is quite possible that the baryon asymmetry, one of the 
most important properties of our Universe, was generated in this transition. Later 
happened the phase transition from the quark-gluon plasma to the hadron matter. 
At the quark-hadron transition the size of a causally connected region of space, the 
horizon radius, was ten kilometers and the age of the Universe 10~^ seconds. 

The new Z(3) phase transition suggested in this work is situated to a temperature 
or energy two orders of magnitude above the electroweak scale. At that time the 
Universe was roughly 10"^^ seconds old and the horizon radius was of the order of 
one micrometer. This hypothetical phase transition is caused purely by the high- 
temperature properties of strong interactions. 

The discussion of the different phase transitions is based on the assumption that 
they are first order. It is pointed out that especially the onset of a cosmological phase 
transition shows a universal behavior. General methods are presented, applicable to 
an analysis of the sequence of events taking place in cosmological first-order phase 
transitions. 

An equation of state is derived for the electroweak matter near the phase transition 
point. The thermodynamically allowed region for the velocities of the phase transition 
front is determined. 

The nucleation rate of bubbles of the broken-symmetry phase is computed in 
generic first-order cosmological phase transitions. The radial dependence of the 
Helmholtz free energy of the bubbles is also discussed. 
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1 Introduction 



Cosmological phase transitions offer a ricli variety of physical phenomena for investi- 
gation, and some of their effects may be observable in the present Universe. In this 
Thesis the mechanisms of cosmological phase transitions are studied, concentrating on 
transitions which are related to strong and electroweak interactions. 

All the phase transitions to be discussed are — or are assumed to be — first order. 
In the first-order transitions a metastable phase may exist alongside a stable one for 
some temperature range. In the thermodynamical limit, i.e., in an infinitely large 
system in which the temperature is changed with an infinitely slow rate, the phase 
transition takes place at a certain critical temperature Tc. 

If the rate of the temperature decrease is finite, as is the case in the expanding 
Universe, the phase transition temperature differs from the equilibrium value Tc. At 
the critical temperature of a phase transition nothing happens, the high-temperature 
phase just moves into a supercooled state. At a somewhat lower temperature bubbles 
of the new phase begin to nucleate. The bubbles grow and convert the space to the 
new phase. The new phase has a lower energy density than the old phase. This means 
that in the phase transition the Universe is heated up to a certain temperature not 
higher than T^. After the transition is completed, the Universe starts to cool again in 
the usual way. 

The supercooling is crucial for the scenarios in which the baryon asymmetry of the 
Universe is generated at the electroweak scale. There is more matter than antimatter 
in the Universe, and more than a billion photons for every baryon. This fundamental 
cosmological fact, crucial for us human beings, should be explained in a satisfactory 
way. 

A few years ago it was realized that at temperatures above the critical tem- 
perature of the electroweak theory, certain electroweak processes mediated by the 
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so-called sphalerons destroy any pre-existing baryon plus lepton number asymmetry 
[Kuzmin, Rubakov and Shaposhnikov 1985]. It became thus necessary to understand 
how the baryon asymmetry of the Universe could have been created at the electroweak 
phase transition. In principle, this is possible since all the three necessary conditions 
for generation of baryon asymmetry [Sakharov 1967] could have been satisfied during 
the transition. Firstly, due to the anomaly in the electroweak theory ['t Hooft 1976a, 
1976b], baryon- number violating reactions were taking place. Secondly, CP-symmetry 
was violated because of fundamental gauge and Higgs interactions of quarks. The 
third condition, a departure from thermal equilibrium, was well satisfied because of 
the supercooling, provided the phase transition is first order. 

In order to obtain any quantitative estimates for the amount of baryon asymmetry 
created, one must have a detailed understanding of how the electroweak phase transi- 
tion proceeded. Motivated by this, we have investigated in this Thesis mechanisms of 
the electroweak phase transition. 

The other phase transitions considered in this work are related to quantum chro- 
modynamics. The possible observable consequences of the cosmological quark-hadron 
phase transition are due to the density inhomogeneities produced during the transition. 
If the length scale of these inhomogeneities had been large enough, they could have 
later affected the nucleosynthesis. This effect could be observed in the abundance of 
light elements in the present-day Universe. However, it seems probable that the length 
scale was too small for that, as will be discussed later on. In the case of the Z(3) 
phase transition suggested in this Thesis, it might in principle be that the density in- 
homogeneities generated could have affected later processes like the electroweak phase 
transition. 

In addition to those questions considered in this work there are several other 
interesting topics related to cosmological phase transitions. For example, we have 
not studied the possible phase transition of a grand unified theory, which may have 
been cosmologically important as a driving source of infiation. Likewise, the topological 
defects, like monopolcs or cosmic strings, which could have been created in cosmological 
first-order phase transitions are not discussed. 
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This Thesis is organized as follows. We first give a brief summary of the contents 
of the original research papers, which are appended. In Section |^ the main events in 
the evolution of the Universe are described. In Section ^ cosmological first-order phase 
transitions are discussed on a general level, without specifying the physical model. In 
Section ^ the general methods presented in the previous Section are applied to two 
physical cases, to the electroweak and the quark-hadron phase transition. Finally, in 
Section ^ we present conclusions and point out directions for the future work. 



Summary of the Original Papers 

Paper I: Nucleation and Bubble Growth in a First-Order Cosmological 
Electroweak Phase Transition. In this paper the thermodynamical properties of 
electroweak matter near the critical temperature are systematically investigated for 
the first time. Assuming a quartic form for the Higgs potential (to be discussed in 



Subsection |4.1| of this introductory review part), we derive an equation of state that 
describes the electroweak phase transition, and compare the electroweak transition with 
the quark-hadron transition. The nucleation rate of bubbles of the broken-symmetry 
phase is computed by solving numerically the field equation. We present a useful ex- 
pression for the volume fraction not touched by the bubbles, slightly different from 
those given previously by other authors. We perform numerical simulations of bub- 
ble nucleation and growth which confirm our analytical calculations. Finally, we also 
study what velocities of the phase front are allowed assuming only that the general 
conditions of energy-momentum conservation and entropy increase are valid. Based on 
these considerations, we claim that in the cosmological electroweak phase transition 
the bubbles most likely grew as weak deflagrations. 



Paper II: Cosmological QCD Z(3) Phase Transition in the 10 TeV Temper- 
ature Range? Using as a starting point the earlier observation that in the QCD 
there are metastable vacua at high temperatures, we develop a cosmological scenario 
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which leads to a phase transition, not known before, at a temperature two orders 
of magnitude above the electroweak scale. Qualitatively, this phase transition differs 
from the usual ones in that the pressure difference between the stable and metastable 
vacua is huge, and in that there were only relatively few bubbles nucleated inside the 
horizon. Our scenario is based on the hypothesis that at very early times domains of 
metastable vacua were created and underwent an inflationary expansion due to some 
processes which could be related for instance to the breaking of the grand unified sym- 
metry. This hypothesis is the main uncertainty in our scenario. Later on it has been 
also claimed that the metastable vacua should only be interpreted as field configura- 
tions that contribute to the Euclidean path integral, not as physically accessible states 
[Belyaev et al. 1992; Chen, Dobroliubov and Semenoff 1992]. In a more recent investi- 
gation it has been, however, argued that the metastable vacua do represent physically 
realizable systems [Gocksch and Pisarski 1993]. 

Paper III: Bubble Free Energy in Cosmological Phase Transitions. In this 
paper the free energy of spherical bubbles is studied in order parameter or Higgs field 
models having the same quartic potential as used in Paper I. A numerical function 
with a good accuracy is given for the nucleation action. Using this nucleation free 
energy of critical bubbles as an input, the general free energy is solved as a function of 
the bubble radius and the temperature. The calculation is based on the approximation 
that all the temperature dependence in the free energy comes from the volume term. 
This approximation should be valid if one is not too far from the limit of small relative 
supercooling. The bubble radius and curvature-dependent interface tension are dis- 
cussed in detail. The results of this study are apphcable for the case of the electroweak 
phase transition, and probably for the quark-hadron transition as well. 
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2 Short History of the Universe 



The big bang model provides a general framework for describing the evolution of 
the Universe. The success of the model is based on only a few — but fundamental — 
astronomical observations: the redshift in the spectra of distant galaxies, the existence 
of the 2.7 Kelvin cosmic microwave background radiation, and the abundance of light 
elements in the Universe. The redshift is believed to be caused by the expansion 
of the Universe. Additional evidence for the big bang model was presented in April 
1992, when it was announced that the observations from the COBE satellite show an 
anisotropy of AT/T fti 6 x 10~^ in the background radiation [Smoot et al. 1992]. 

In principle, we are able to obtain (semi) direct information of the early Universe 
by observing electromagnetic or gravitational radiation, or exotic relics like very mas- 
sive particles or small black holes. So far, the only cosmological messengers we have 
been able to detect in our instruments are the photons. This means that in direct ob- 
servations we have to limit ourselves to a Universe older than a few hundred thousand 
years. At that time the temperature was somewhat less than the binding energies of 
electrons in light atoms. The electrons and light nuclei were able to form stable atoms, 
and the Universe became transparent for photons. This is the epoch we are looking at 
when observing the cosmic microwave background radiation. 

The abundance of the light elements, as observed in the present-day Universe, 
gives us indirect but firm evidence of the time of the primordial nucleosynthesis. The 
major part of the nucleosynthesis took place when the Universe was a few minutes old 
[see e.g. Applegate, Hogan and Scherrer 1987]. This is the earhest epoch of which we 
have more or less certain information. 

The study of the very young Universe requires an extrapolation of the cosmological 
model to even earlier times. In addition to Einstein's theory of gravitation, what is 
needed for that is a knowledge of interactions between elementary particles at high 
energies. This knowledge, the standard model of particle physics, is based on laboratory 
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Figure 1: Evolution of the Universe. 



experiments done at colossal particle accelerators. 

It is believed that when the Universe was approximately 10~^ seconds old and 
temperature was of the order of 100 MeV, a phase transition from the quark-gluon 
plasma to the hadron matter took place. In the quark-gluon plasma phase the quarks 
and gluons were free, whereas in the hadron phase they became confined and formed 
mesons like the pions, and baryons like the proton and the neutron. The existence 
of quark-gluon plasma has not yet been confirmed experimentally, but several groups 
are trying to detect it by using heavy-ion collisions. As mentioned in the Introduc- 
tion, probably the scale of the density inhomogeneities produced in the cosmological 
quark-hadron phase transition was too short to have observable consequences (see 
Subsection [4.21 ). 

Further back in time, the electroweak phase transition took place at t ~ 10^^^ s. 
The temperature was then about 100 GeV which roughly corresponds the highest 
energies available in particle accelerators. Hence this transition presents the earliest 
time in the history of the Universe of which we have a moderately detailed knowledge. 
At the classical level of the electroweak theory, the mass terms of the quarks and 
gauge bosons vanish in the high-temperature symmetric phase. In the phase transition 
their mass terms, which are proportional to the vacuum expectation value of the Higgs 
field, become non-vanishing. It is believed that the electroweak phase transition had a 
significant effect on the baryon asymmetry observed in the present Universe. 
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Of still earlier times we presumably know at least how temperature evolved during 
most of the time. Extrapolating the theories of strong and electroweak interactions up 
over ten decades in energy, one is tempted to believe in grand unification at the scale 
of 10^^ GeV. However, the huge energy gap between the electroweak and the grand 
unified scales may very well hide new phenomena. One of the most speculated issues is 
supersymmetry; if it exists, the symmetry between bosons and fermions was restored 
above certain temperature. The hypothetical transition suggested in Paper II is also 
situated to the gap between the electroweak and grand unified scales; its transition 
temperature is two orders of magnitude above that of the electroweak phase transition. 

The most important idea, connected to the (hot) big bang model and related to 
these very early times, is the hypothesis of inflation. The inflation, or an exponen- 
tial expansion of the scale factor at very early times, would resolve the smoothness 
problem — why the temperature of the microwave background is almost uniform over 
scales bigger than the horizon scale when the photons last scattered. Furthermore, the 
inflation predicts that the Universe expands eternally but with an always slowing rate, 
the alternative we are astonishingly close to according to observations of the average 
mass density of the present Universe. 

The extreme limit for the validity of the standard cosmology is the Planck scale. 
Quantum mechanics tells that if we inspect any system on a scale small enough, the 
physical quantities fluctuate strongly. For gravity this fluctuation length scale is given 
by 1/Mpi, where the value of the Planck mass is Mp\ ^ 1.2 x 10^^ GeV.0 The horizon 
radius, the size of a causally connected region of space, was equal to this fluctuation 
scale at the Planck time when the age of the Universe was 10-^"^ s. At the Planck scale 
general relativity is not valid any more, one would need a quantized theory of gravity. 

Let us finally present the fundamental equations governing the evolution of the 
Universe in the big bang model. A more detailed treatment can be found for instance 
in the textbook by Kolb and Turner [1990]. A beautiful way to formulate the Einstein 

""^In this Thesis the system of so-called natural units is used: fi — c^ks = However, numerical 
values for time and length are often given in Sl-units. 
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equations is to use the action 

J d'^x^/^ {^A,<^,v, - Yg^^sc} , (1) 

where g is the determinant of the metric tensor, Ca,4>,'^ is the Lagrangian density for all 
the gauge and matter fields in the standard model, G is the gravitational constant and 
Rsc the curvature scalar. The Einstein equations now follow by demanding that the 
variation of the action with respect to the metric vanishes. The metric that produces 
homogeneous and isotropic three-spaces is the Robertson-Walker metric, 

df^ 



+ f^d6'^ + sin"" e d</7 



1 — /cf" 



(2) 



where R{t) is the cosmic scale factor, f is a dimensionless scaled coordinate and the 
constant k takes values of +1, or —1 for three-spaces of positive, zero or negative 
curvature, respectively. 

The full stress-energy tensor, the variation of the action of CA,(j>,ip with respect 
to the metric, is a complicated object. However, to obtain the general evolution of 
the Universe it is enough to assume that it is given by the simple perfect fluid form 
T'^y = diag(£, — p, — p, — p), where e is the energy density and p the pressure. The 
independent Einstein equations can in this case be chosen as 

d{sR^) = -pd{R^). (4) 

The upper one, or Priedmann equation, gives the time dependence of the scale fac- 
tor. The lower equation can also be derived from the local conservation of energy- 
momentum of the perfect fluid. 

In the case of the very early Universe, when matter behaves like radiation, pressure 
and energy density are given by 

^ 3 90^ ^ ^ 

where is the effective number of relativistic degrees of freedom. Every relativistic 
bosonic and fermionic degree of freedom increases the value of g^, with unity and 7/8, 
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respectively. The relation (^) is an idealization, valid for free particles. Inclusion of the 
interactions between particles modifies it slightly even at temperatures far away from 
any phase transition [see e.g. Enqvist and Sirkka 1993]. 

When the early Universe was at local thermal equilibrium, entropy was conserved. 
This can be seen from eq. (0) using standard thermodynamics combined with the 
extreme smallness of all chemical potentials. The conservation of entropy implies that 
RT = constant, as long as g^, does not change. The curvature term (the one with k) 
is negligible in the Friedmann equation (^) and we can neglect it. Now it follows from 
substituting the expression for the energy density in eq. (H) to the Friedmann equation 
and expressing the gravitational constant as G = 1/M|i that 



This important equation gives the relation between temperature and time in a 
radiation-dominated Universe. 




(6) 
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3 First-Order Phase Transitions 



3.1 Nucleation Rate in Field Theory 

Quantum field theory forms tlie tlieoretical framework for describing microscopic rela- 
tivistic processes at zero temperature. Tlie theory is elegantly formulated in terms of 
the Feynman path integrals. This formulation can be generalized to non-zero temper- 
atures [Bernard 1974; Gross, Pisarski and Yaffe 1981]. 

The thermodynamical properties of a physical system at nonzero temperature can 
be calculated from the partition function 

Z = tT e-^^ , (7) 

where f3 is the inverse temperature and H the Hamilton operator. 

The operator exp(— in eq. (0) is analogous to the quantum- mechanical time 
evolution operator exp(— iift), if one makes the identification it = 1/T. Indeed, in 
Euclidean space-time the path-integral representation for the partition function of a 
scalar field is 

Z = f [d0(r,x)] expj- /%r / d^x £^^(0, 90)1 , (8) 

J/3-periodic [ JO J J 

where the imaginary time is denoted by r and Ce is the Euclidean Lagrangian density. 
We can see that along the imaginary-time direction the field propagates over a 
distance (3. Because of the trace in the definition of the partition function in eq. 
the physical states have to be identical at "times" and (3. Thus the scalar field must 
obey 

0(r + /5,x) = 0(r,x) Vr,x. (9) 

The subscript periodic' in eq. (H) refers to this periodic boundary condition. 

For fermion fields the boundary conditions are not periodic but antiperiodic. This 
is due to the anticommuting nature of fermions, which is reflected in the definition of 
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the time-ordering operator. For gauge fields Gauss's law (known in electrodynamics as 
Maxwell's first law) must be imposed on physical states. Furthermore, periodicity of the 
temporal gauge field Aq is just the most natural choice, not an automatic consequence. 

At this point, two comments related to the partition function are in order. Firstly, 
by formulating the theory in Euclidean space-time we loose information on the real-time 
evolution of the system. The imaginary-time formalism is applicable only for describing 
processes which occur at or near thermal equilibrium. Most of the time the evolution 
of the Universe probably happened very near thermal equilibrium, but not during the 
first-order phase transitions. However, the moment of the onset of a first-order phase 
transition can still be calculated using these methods. Secondly, we have assumed that 
the chemical potentials for different particle species vanish. This is usually justified 
when one calculates quantities like pressure or energy density because of the very small 
baryon and (presumably) lepton asymmetry [see e.g. Kajantie and Kurki-Suonio 1986]. 

Next, we will discuss the nucleation rate in phase transitions. We consider a single 
scalar field, for which the Euclidean effective action is given by 



The potential V{(f), T) should be understood as an effective potential. Besides the classi- 
cal part, it includes both zero-temperature quantum corrections and finite-temperature 
thermal corrections. We assume that the corrections do not change the derivative part 
of the action. For a discussion on this point see for instance [Brahm 1992] . 

In fig. 2, the qualitative behavior of the potential is shown. At high temperatures, 
the only vacuum of the scalar field model is at = 0. At some lower temperature, 
another local minimum starts to develop. With decreasing temperature, the local min- 
imum soon becomes the true vacuum, and the vacuum at origin becomes metastable. 
At still lower temperatures the metastable vacuum vanishes, and the system has again 
only one vacuum. It is the wall between the old metastable vacuum and the new 
true vacuum that causes the phase transition to be first order. The thermodynamical 
transition temperature T^. is the temperature at which the two vacua arc degenerate. 
In cosmology the system supercools because of the expansion of the Universe and be- 
cause the transition rate is only finite, and therefore the transition takes place at a 




(10) 
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Figure 2: Qualitative behavior of the effective potential for the scalar field, 
F(0,r). The potential is plotted at four different temperatures which obey 
T4 < T3 < T2 < Ti and is normalized to vanish at ^ = for all values of T. 

temperature somewhat lower than Tc- 

We assume that the mechanism for the phase transition is homogeneous nucle- 
ation. However, one should note that the presence of relic fluctuations or exotic objects 
like magnetic monopoles, cosmic strings, black holes or very massive particles might 
modify the mechanism. In the phase transition, bubbles of the stable phase arc cre- 
ated. Nucleated bubbles which are bigger than so-called critical bubbles begin to grow. 
When the bubbles grow, the Universe is gradually converted to the new phase. 

The decay rate of a metastable vacuum has in the context of statistical field theory 
been calculated by Langer [1969]. For relativistic quantum field theory at zero temper- 
ature the decay rate has been evaluated by Coleman [1977] and Callan and Coleman 
[1977], and the same semiclassical methods can be applied also at finite temperatures 
[Affleck 1981; Linde 1977, 1981, 1983; Arnold and McLerran 1987]. At high tempera- 
tures, the mechanism for bubble creation is, instead of quantum tunnelling, thermal 
over-barrier nucleation. The four-dimensional action can be approximated with the 
three-dimensional one: 



In the case of weakly first-order phase transitions the three-dimensional action is at 
least a very good approximation, if not even exact. 




(11) 
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The so-called critical bubble is a non-vanishing solution of the field equation: 

=-V^0 + mT) = O, (12) 

where the prime means derivative with respect to 0. The boundary conditions are 
that the derivative of the solution must vanish at the center of the bubble, cho- 
sen to be at the origin, and that at infinity the solution must be in the metastable 
vacuum. It has been shown that the critical bubble is spherically symmetric 
[Coleman, Glaser and Martin 1978]. An example of the critical bubble solution 0(r) is 



given later on in fig. 6 in Subsection 4.1 



The decay rate of the metastable vacuum can be calculated from the imaginary 
part of the partition function. The extremum action corresponding to the critical 
bubble will be denoted by 6*3 (T), and the potential is normalized as in fig. 2 so that the 
vacuum action of the high-temperature phase vanishes. The probability of nucleation 
per unit time per unit volume is 



piT) = — 



lS3{T)Y^^ det'[-d^ + V"{^{r),T)] 



2n \ 2nT J 



exp 



S,{T)/T , (13) 



det[-d^ + V"{0,T)] 

where det' means that the zero eigenvalues resulting from the three translations of the 
bubble center should be omitted when calculating the functional determinant. The 
real quantity U- is the (angular) frequency of the unstable mode of small fluctuations 
around the bubble solution. 

The dominant factors in the nucleation rate (0) are the exponential part and the 
dimensional part of the prefactor. By prefactor we mean the product of all the factors 
that multiply the exponential part. If one transforms the functional determinant to di- 
mensionless units, the dimensional quantity that factorizes out is M^(T), where M(T) 
denotes the mass of the quadratic term of the potential in the symmetric phase. It has 
been shown [Brihaye and Kunz 1993] that the frequency of the unstable mode can be 
approximated well — at least for some set of parameters of the quartic potential — with 
the thin-wall formula cu^ = -\/2/-RTcr, where -Rtct is the radius of maximal tension of the 
critical bubble (see Paper III). The remaining dimensionless determinant is estimated 
to differ from unity at most by a few orders of magnitude. This expectation gains some 
confidence from a somewhat similar case, namely from investigations of the sphaleron 
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transition rate where the corresponding determinant has been calculated both numeri- 
cally [Carson and McLerran 1990] and analytically [Carson et al. 1990]. Very recently, 
the fluctuation determinant has been evaluated in the case of a critical bubble as well 
[Baacke and Kiselev 1993]. The result, though renormalization scheme dependent, is 
not very far from unity for realistic values of T/M. 

Changing the value of the prefactor with even several orders of magnitude would 
not signiflcantly affect the cosmological nucleation calculations, as we shall see later 
on. We may rewrite the nucleation rate in the form 

where the value of the slowly changing function h{T) is not essential compared with 
the exponential part. Besides the simple estimates and T^, prefactors like M^(T) 
[McLerran et al. 1991] and the Laplace radius of the critical bubble i?Lcr multiplied 
by some other dimensional factors [Csernai and Kapusta 1992a] have been proposed. 
(i?Lcr will be deflned in eq. (|30| ) in Subsection |373| .) 



3.2 Bubble Nucleation and Growth 

A measure telling how the phase transition proceeds is the fraction of space con- 
verted to the new phase. In cosmology, the formula for the volume fraction remain- 
ing in the old phase given by Guth and Tye [1980] is usually employed [see also 
Guth and Weinberg 1981]. For time scales much shorter than the Hubble time the 
expansion of the Universe can be neglected. This approximation should be valid for 
the whole electroweak phase transition and for the initial stages of the quark-hadron 
phase transition. The fraction of space still in the metastable phase at time t is given 

by 

Uit) = e-^W; J(t) = fdt'piTit')) V{t',t) , (15) 

where tc is the time when temperature is equal to Tc and V{t', t) is the volume that a 
bubble nucleated at time t' occupies at time t. 
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A different expression for /ms has been proposed by Csernai and Kapusta [1992b]. 
In their approach, the metastable fraction is given by the integral equation 

Us{t) = 1 - fdt'p{T{t')) U{t') V{t',t) . (16) 

The two expressions (|15]) and ([T6|) coincide in the very beginning of the phase 
transition when 1— /ms ^ 1? and always if the bubbles do not grow. However, they give 
different results when growing bubbles begin to overlap. The underlying assumption 
of the Guth-Tye formula (p!5| ) is that when two bubbles come into contact they simply 
grow inside each other without any effect on those parts of the bubbles which do 
not touch yet. As a consequence the phase transition is never fully completed in an 
infinite volume; there always remains a non-zero fraction of space in the old phase. The 
exponentiation of the naive fraction I{t) takes care of bubble overlap and of the fact 
that the volume fraction where new bubbles can nucleate decreases with time. The 
latter feature is explicitly taken into account in the Csernai-Kapusta formula (p!6[) , 
whereas the volume occupied by two bubbles which have come into contact with each 
other is still given as a sum of their individual volumes. 

Before going further, the concept of bubble in this context should be clarified. 
Nucleated bubbles grew probably as deflagrations both in the electroweak and in the 
quark-hadron phase transition. The effect a deflagration bubble has on the cosmic 
fluid is illustrated in fig. 3. The coordinate system is chosen to be initially at rest with 
respect to the cosmic fluid. First, the fluid is hit by a supersonic shock front which 
heats it and accelerates it to the velocity f fluid in the same direction as the shock. Later 
on, a subsonic deflagration front, which is the actual phase boundary, propagates into 
the moving fluid transforming it to the new phase. The fluid is cooled down and stays 
again at rest. Depending on the situation, the bubble surface should be identified in 
the calculation as either the phase boundary or as the shock front. In the latter case 
the expressions (p!5| ) and (p!6| ) give, instead of the fraction of space in the old phase, 
the fraction not touched by the shocks. 

When bubbles begin to touch several things can happen. If the deflagration 
fronts expand very slowly, as is the case in the quark-hadron phase transition after the 
system has reheated back to Tc (Subsection |4^ ), they rearrange the surface quickly 
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Figure 3: World lines of cosmic fluid near a deflagration preceded by a shock 
front. In this idealized flgure space-time is 1+1 dimensional, and the initial 
bubble, located at the origin, is infinitesimally small. Dashed lines show how two 
fluid points move. In the old phase, temperature is Tf outside the shock, and Tq 
in the region between the shock and the interface. In the new phase temperature 
is equal to Th. Slopes of the lines are the inverse velocities of the shock front, the 
phase boundary ("wall"), and the fluid in the intermediate region. 

to a spherical shape so that the total volume stays approximately constant. However, 
when the radii of the bubbles exceed a characteristic length Rf^^s{t), the rearrangement 
process is too slow compared with the cosmic expansion to take place [Witten 1984]. 
Finally, when the stable phase fills most of the space, the interface has arranged itself 
in such a way that the high-temperature phase is located in isolated droplets. For slow 
deflagration fronts in the electroweak phase transition the process should be similar, 
except that the characteristic length Rfus{t) is smaller than the one which follows 
directly from Witten's argument since the time scale of the whole phase transition is 
much shorter than the Hubble time (Subsection [4.1| ). The case of shock collisions which 
is relevant for studying the reheating is qualitatively different, because shock fronts do 
not fuse together. 

When comparing the Guth-Tye and Csernai-Kapusta expressions for the volume 
fraction fms{t), one may come to the conclusion that neither of them is perfect. For 
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phase boundaries, the Guth-Tye formula overestimates the time needed to complete the 
phase transition because it assumes that the bubbles grow inside each other without any 
interference; the Csernai-Kapusta formula underestimates the completion time since 
after two bubbles have fused together the single volume grows in reality more slowly 
than the total volume of two bubbles had they been separated. For shock fronts, the 
Guth-Tye formula should be rather accurate unless reflection plays a significant role 
in the front collisions. In actual calculations the Guth-Tye formula is easier to handle. 
We will adopt the Guth-Tye formula (|15D for further use. 

Let us now inspect more closely that stage of the phase transition during which 
most of the bubbles are nucleated. We will assume that the nucleation is not active 
any more inside the bubbles. This assumption is too strong only if we consider slow 
deflagrations preceded by shocks, and if the shocks are, on one hand, too weak to reheat 
the plasma enough to stop the nucleation and, on the other hand, not so weak that 
their effect on the nucleation rate could be neglected. In the last case, if one identified 
the bubble boundary as the defiagration front and not as the shock, there would be no 
complications. 

The actual period of nucleation is short. During it the prefactor in the expres- 
sion ([l^) for the nucleation rate changes only slowly compared with the exponentiated 
nucleation action, as does the bubble wall velocity v. We may hence approximate them 
as being constants. By writing the nucleation rate in the form 



The size of the bubble when it is first nucleated is very small, as will be shown later on, 
and therefore it is left out of this expression. We define t/ as the time when the bubbles 
occupy a significant fraction of the space. For simplicity tf is always called the phase 
transition time, even if the bubbles are defined in terms of the shock spheres. The above 
expression for /ms can be evaluated by noticing that all the essential change in it takes 
place when {t — tc)/{tf — tc) is close to unity [Fuller, Mathews and Alcock 1988]. After 



p{t) =poe 



Sit) 



(17) 



we obtain the following expression for the metastable volume fraction: 




(18) 
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expanding the nucleation action as a Taylor series, the naive fraction I{t) which after 
exponentiation gives the fraction of space still in the old vacuum, exp[— /(t)] = /ms(^), 
can be expressed as 

I{t) = I{tf) e-^(*/-*) , (19) 

where 

The positive factor (3 turns out to be an important scale-setting parameter for the 
phase transition. 

The cosmological phase transition discussed in Paper II differs qualitatively from 
the 'normal' cases. It is exceptional in the sense that in thermal equilibrium the same 
phase is the true vacuum at all temperatures. The approximation presented in last 
paragraph is valid both for this 'exceptional' phase transition and for the electroweak 
and the quark-hadron transition when 

/3(t-g>i. (21) 

In normal transitions, the nucleation action is proportional to the square of l/{t—tc), 
and decreases very rapidly after the critical temperature has been reached. As a 
consequence, the validity condition is very clearly fulfilled in both electroweak and 
quark-hadron transition: the left-hand side of eq. (|2T| ) is during actual nucleation 
bigger than unity by two orders of magnitude. In the transition considered in Paper II, 
the nucleation action behaves as log^''^(l/t). (In this case tc in eq. (pTj) means the time 
when the metastable vacua were created.) Due to the slowness of the decrease of the 



action, the left-hand side of eq. (^) is not more than about 6. However, this value can 
still be considered to lie within the allowed range. 

As in Paper I, we will define to be the time when the volume fraction occupied 
by the bubbles equals 1/e, in other words, we put /(t/) = 1 in eq. (plQ]). The phase 
transition time tf can then be solved from the equation 

1 = m) = ^pitf) ■ (22) 

In principle, the Guth-Tye formula may not be any more fully valid at t/ due to bubble 
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Figure 4: A two-dimensional simulation of bubble growth, taken from Pa- 
per I. The three frames of size (40f//3)^ show the bubble configuration at times 
tf — 4.5//3, tf — 2.5/13 and tf — 0.5/13. Bubble collisions were neglected in this 
simulation, but new bubbles were not allowed to nucleate inside existing ones. 



collisions. However, the inaccuracy that this causes is negligible when determining for 
example how much the system supercools. 

As an interlude, let us approximate eq. (^) dimensionally as p{tf)t\ 1. The 
purpose of this zeroth-order approximation is to give an estimate of the scales involved. 
Utilizing the relation between time and temperature in radiation-dominated Universe, 
eq. (0), we obtain 



5(t/) = ^=^^41og^, (23) 

where Tf stands for the temperature in those parts of space which have not been affected 
by the bubbles yet (see fig. 3). Now we can see why uncertainties of several orders of 
magnitude, e.g. in the prefactor of the nucleation rate, are not important in cosmology: 
the value of the critical nucleation action is very large — as a result of the very slow 
expansion rate of the Universe or the weakness of the gravitational interaction, which 
is seen as the large value of the Planck mass. For the transition temperatures of the 
electroweak and the quark-hadron phase transition, the right-hand side of eq. (p3D is 
equal to 150-160 or 180-185, respectively. This also means that the validity of the 
semiclassical or WKB approximation, employed in the calculation of the nucleation 
rate, should in this regard be on firm footing. 
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Let us return to the general, more accurate analysis. As in Paper I, we define an 
effective nucleation rate ip{t) as follows: 

m=Pit)fut)--'-f'-^. (24) 

The decay rate per unit volume per unit time is corrected with the metastable fraction, 
since as the phase transition proceeds there is less space available for the bubbles of the 
new phase to nucleate. Integrating the effective nucleation rate we obtain the number 
of bubbles in unit volume function of time: 

nit) = I' dt'm = ^ [1 - /-(^)] • (25) 

The number density increases linearly with the fraction of space in the low-temperature 
vacuum. 

A good estimate for the final number density of bubbles nfinai is obtained by setting 
/ms = in eq. (|25|) and using eq. (^). If several collisions of the shocks are needed to 
produce notable reheating (Subsection |4.2| ), the true final number density is somewhat 
larger than the one obtained from eq. ([25|). The average distance of nucleation centers 
Rnnch defined as nfinai = l/-Rnucb is given by 

i?nuci = 27ri/3^ , (26) 

where (3 is defined in eq. (^Of). Increase in the rate of change of the action means a 
decrease in the nucleation distance because then more bubbles nucleate during a given 
time interval after the nucleation has effectively been turned on. 

In both the electroweak and the quark-hadron phase transition there were a vast 
number of bubbles nucleated inside the horizon. This was partially due to the rapid 
change in the nucleation action, and partially due to the smallness of the dimensionless 
ratio ct^/^/La/T^, where a is the interface tension and L the latent heat of the transition 



(see eq. ( pT]) in Subsection p.3[ ). On the other hand, in the case of the phase transition 
suggested in Paper II the nucleation action decreases so slowly that there were only a 
few hundred nucleated bubbles inside one horizon. 

Let us now consider the distribution of bubble sizes. This has been recently dis- 
cussed by Turner, Weinberg and Widrow [1992]. The bubble size distribution is easily 
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Figure 5: Bubble size distribution. The area to the right of the right vertical line 
is the bubble size distribution at time t = tf — 3/P, and the area to the right of 
the left vertical line is the bubble size distribution at time tf. In other words, 
the vertical lines correspond bubbles of zero-radius at those times. 



obtained after realizing that a bubble which has a radius p at time t, was nucleated 
at time t—p/v. Here we made again the approximation that at nucleation the bubble 



radius is so small that it can be neglected (see Subsection p.3|) . The distribution of 
bubble sizes is hence given by 

g,{p) = ^^P(t-p/v). (27) 

V 

The number of bubbles in unit volume at time t with radius between p and p + dp is 
equal to gt{p)dp. More explicitly, the bubble size distribution is given by 

gtip) = ^ e-" exp [-e""] ; u = /3itf-t + p/v) . (28) 

When inspecting the spatial scale of density inhomogeneities, the distribution of dis- 
tances between nucleation sites of bubbles near each other could likewise be used. 
Meyer et al. [1991] have studied it extending the naive fraction approximation for the 
whole nucleation period. 

The distribution of bubble sizes given in eq. (^) is presented in fig. 5. When time 
t is in the vicinity of the phase transition time t/, our approach loses its reliability. 
Therefore the distribution at time t/, the area to the right of the left vertical line in the 
figure, may well be somewhat distorted near the origin, whereas the other distribution 
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corresponding to t = t/ — (3//3) can be expected to be quite accurate. The bubble 
size distribution ( P^ ) is universal in the sense that if the validity of the assumptions is 
equally good in two different phase transitions, the appropriately scaled distributions 
should look the same. 



3.3 Thin- Wall Limit 

In the so-called thin- wall limit, or small relative supercooling limit, the equations for 
the amount of supercooling and the distance between nucleation sites can be expressed 
in a simple form and solved with a good accuracy. Let us consider a thin-walled bubble 
of radius R. Let the whole interior of the bubble be in the vacuum of the new phase and 
the whole exterior in the old vacuum, and let the interface between be infinitely thin. 
The free energy density difference between the two phases is, in the absence of any 
relevant conserved charge, equal to minus the pressure difference Ap. By difference 
we mean the value of the quantity in the low-temperature phase minus that in the 
high-temperature phase. Hence the free energy of the bubble is 

F{R) = -—ApR^ + AtxgR^ , (29) 

where a is the interface tension. A generalization of this expansion is discussed in 
detail in Paper III. The radius of the critical bubble and the corresponding free energy 
density are found by maximizing F{R) with respect to i?, with the result 

Above i?Lcr is Laplace's definition for the radius of the critical bubble. The probability 
of fluctuation in which a critical bubble is formed is proportional to the Boltzmann 
factor exp(— Fcr/T). 

For the normal phase transitions the pressure difference can be written as Ap{T) ^ 
L(l— T/Tc), as given in Section ^ in eq. (P7|). Using eq. the nucleation action 



appearing in eq. ([T7|) can be expressed as follows: 
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where L is the latent heat of the transition. In the case of the exceptional phase 
transition discussed in Paper II, one had to use the exact expression for the pressure 
difference. Thanks to the simple form of the expression this would not give rise to any 
difficulties in the thin-wall calculation. In the rest of this Subsection, we will analyze 
only the normal phase transitions. 

Utilizing eq. @ and using the fact that T is very close to T^, we can write the 
derivative (3 of the action in eq. (|20|) as tc(3 = S(Tf)/{l —Tf /Tc). Substituting the 
expression (|1^ for the nucleation rate, the rightmost part of eq. ( p2D then becomes 



where 



■ e " // = exp 



Tf/T,Y ^ (1 - ry/r,)3/2 



(32) 



A ^ 41og (^) + 21og (^^) + log(8..3,) , (33) 

and h is the dimensionless factor introduced in eq. (p!4D. The phase transition tem- 
perature is obtained by demanding that the argument of the exponential function on 
the right-hand side of eq. (|3^ ) vanishes. A good accuracy is achieved by iterating the 
resulting equation once. This procedure gives the following expression for the amount 
of supercooling [Kajantie 1992]: 

1-^ = ^; i = A- 4 \og{A'/'/C) . (34) 

J-c va 

The phase transition temperature Tj given by this equation should be close to the 
equilibrium transition temperature Tc, otherwise the approximation used is invalid. 

In the expression for A in eq. (|33| ) the first term is by far the dominant one, while 
for A the second term does have some significance; for example in realistic estimates 
for the quark-hadron transition its value is 30% of that of the first term in the ex- 
pression (plf). The physical reason for such a large value of this correction term is 
that the growth time of bubbles, although tremendously longer than the time scale of 
microscopic interactions, is still very far from the Hubble time. 

Now we can determine the average distance between nucleation centers. We com- 
pare it with the horizon radius i?hor = 2tc: 

6.00 ^- = . (35) 
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The later a phase transition takes place, the more bubbles there tends to be inside the 
horizon. For example, in a phase transition at the grand unified scale there should be 
less bubbles created than in a phase transition at the QCD scale. Of course, in actual 



dependence on the factor log(Mpi/T) related to the cosmic background. 

It is also interesting to calculate the radius of the critical bubble Rcr, i-e., the 
bubble size immediately after the nucleation. It can be determined from Laplace's 
relation (^) and is given by 



Here is compared with l/T^, which gives the length scale for microscopic processes. 

From eqs. ( PBD and ( |5BD we can observe that a natural length scale for -Rnuci is 
roughly set by the horizon radius, and for by the microscopic length 1/Tc. We can 
thus conclude that the bubbles are indeed microscopic when first nucleated, but grow 
to a macroscopic size before colliding with other bubbles (unless the interface tension 
is extremely small). 



cases the dependence of -Rnuci/-Rhor on the microscopic ratio a^^'^/L may well cancel the 




3 A 



(36) 
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4 Phenomenology of Cosmological Phase 

Transitions 



In this Section the general methods described above are apphed to the case of the 
electroweak and the quark-hadron phase transition. We present a partially quantitative 
description of the different events that took place during these cosmological transitions. 

We consider in this introductory review part only the first-order electroweak and 

quark-hadron phase transitions and their cosmological consequences, but whether these 
transitions in Nature really arc first order is by no means known with certainty. A 
second-order phase transition would have had less cosmological significance. 

A fundamental ingredient for the study of a first-order phase transition is the 
thermodynamical equation of state. The absence of chemical potential simplifies the 
thermodynamics, since the physical quantities become functions of only one variable, 
the temperature. The pressure p{T), now equal to minus the free energy density, is 
taken as the basic quantity. 

When the temperature is in the vicinity of the thermodynamical transition tem- 
perature Tc, the pressure difference between the stable and the metastable phases can 
be expanded as 



The quantity Ap(T) is defined as the pressure in the low-temperature phase minus the 
pressure in the high-temperature phase. However, this expansion is not valid for the 
exceptional phase transition suggested in Paper II because no critical temperature Tc 
exists there. 

Given the pressure p{T), the entropy and the energy density can be determined 

from 




(37) 



s{T) 



dp 
dT ' 



e{T) ^Ts-p. 



(38) 
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These relations hold in both phases. 



4.1 The Elect roweak Phase Transition 

Although the concept of symmetry restoration at high temperatures has been 
known for long [Kirzhnits 1972; Kirzhnits and Linde 1972; Dolan and Jackiw 1974; 
Weinberg 1974], the electroweak phase transition has remained poorly understood un- 
til recent times. The interest in it was renewed some time after the observation that 
nonperturbative processes, mediated by the sphalerons, have a significant effect on 
the baryon asymmetry of the Universe [Kuzmin, Rubakov and Shaposhnikov 1985]. 
Since then, several authors have studied the dynamics of the first-order electroweak 
phase transition [McLerran et al. 1991; Turok 1992; Paper 1; Anderson and Hall 1992; 
Dine et al. 1992; Liu, McLerran and Turok 1992; Carrington and Kapusta 1993]. 

The minimal standard model of electroweak interactions contains one scalar dou- 
blet. The order parameter of the electroweak phase transition is the field corresponding 
to the low-temperature physical Higgs particle, usually taken as the real part of that 
component of the doublet which acquires a vacuum expectation value. 

The Higgs boson has never been experimentally detected, and guesses for its mass 
cover a wide range. Even its existence may be questioned; for instance, a condensate 
of heavy quarks could effectively act as a Higgs particle [Lindner 1992]. Futhermore, 
it might be that the real electroweak theory contains more scalar fields than just one 
doublet. In the case of two doublets a finite temperature potential of the same form 
as in the case of one doublet can be used as a reasonable approximation for that 
combination of scalar fields which drives the transition [McLerran et al. 1991] . 

The first-order electroweak phase transition is commonly described by the follow- 
ing effective quartic potential [Linde 1983]: 

y T) = ]p{T^ - T^)cf>' - ^aTcf>' + ^A0^ , (39) 

where the order parameter field 4>{x) is a real scalar function. The qualitative (but 
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not quantitative) behavior of this potential is shown in fig. 2. The thermodynamical 
transition temperature is 

Tc = , ^° , , (40) 

V 9A7 

and the lowest temperature where the symmetric vacuum can exist is Tq. 

The potential V{(f), T) could in principle be derived from the full microscopic 
theory. The question how to best determine the effective potential beyond the naive 
one-loop level is currently under an active study. One method is to employ effec- 
tive three-dimensional theories, with [Bunk et al. 1993; Kajantie, Rummukainen and 
Shaposhnikov 1993; Farakos et al. 1993] or without [Shaposhnikov 1993] lattice Monte 
Carlo simulations. These studies seem to indicate that the electroweak phase transition 
would not be as weakly first order as expected on the basis of perturbative analysis. 
(There are arguments telling that the existence of supersymmetry could also make the 
transition more strongly first order [Espinosa, Quiros and Zwirner 1993].) 



The point of view taken here is that the potential in eq. (|39D should be regarded as 
a phenomenological one, valid in the vicinity of T^. The parameters Tq, 7, a and A are to 
be chosen so that the potential quantitatively correctly describes the phase transition. 
In principle, the values of these parameters could be determined by experiment or 
observation. In practice, first-principles calculations, even if not fully satisfactory, are 
employed to obtain an estimate of the relevant scale of the parameters. 

Next, we will discuss the equation of state near the electroweak phase transition. 
Let us inspect two configurations where the field is spatially uniform with a value 
corresponding either to the minimum at the origin or the other minimum of V{(f),T). 
Since the potential vanishes at origin, the difference in the free energy densities between 
these two configurations is equal to the value of the effective potential at the other 
minimum. We denote it as V"(f (T), T) = —B(T), where v{T) is the value of the field at 
the other minimum. This leads to the following expressions for the pressure (Paper I):i 

p^{T) = aT\ pr,{T) = aT^ + B{T) . (41) 

The same factor a is used in both radiative terms, because at Tc the pressures of both 



^For simplicity the subscripts 'q' and 'h', adopted from the quark-hadron transition, are used also 
here to denote the high- and low-temperature phase, respectively. 
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phases must be equal and B(T) vanishes. The term B{T) was attached to the low- 
temperature phase since both B(T) and the low-temperature phase exist only up to a 
certain temperature somewhat above T^. Moreover, one might add a constant — B(0) 
to each of the phases in order to make Ph{T) vanish at zero temperature. However, this 
is unnecessary since the above equation of state is invalid at low temperatures anyway 
(Paper I). 



From the relation (|g) one observes that the scale for the constant a in eq. (^If) 
is apart from the factor 7r^/90 set by the effective number of relativistic degrees of 
freedom. At very high temperatures all the particles of the minimal standard model 
contribute to g^: giving 

7 



^* =2 +3 X 2 + 8 X 2+ 4 +- 



7 ^wcak g $ 



106.75 . (42) 



3x4+3x2+6x3x4 

e,ii,T V q 

At the phase transition all three weak gauge bosons increase their degrees of freedom 
with unity by eating the altogether three would-be Goldstone bosons from the scalar 
doublet. This does not change the value of significantly, because slightly below 
the masses of the gauge bosons and the physical Higgs particle are presumably small 
compared with the temperature. However, the top-quark, the other yet undetected 
particle predicted by the minimal standard model, might decrease the value of g^, 
significantly — the contribution from a single quark to the value of g^. is about 10%. 

At the phase transition, the jump in g^, folded together with the change in the 
interactions causing deviation from the radiative free-gas behavior, produces latent 



heat of the transition. For the equation of state (^) the latent heat is given by 
L = —TcB'{Tc), and its value is hence reflected in the parameters of the potential (|39D. 
One should note that latent heat L includes not only the effect of the Higgs particle 
but also of all the other particles, since in our order parameter model all the other 
fields have been integrated out. 

As we have now specified the potential V{(f), T) we will shortly return to the 
nucleation process. It is convenient to express the temperature dependence of various 
quantities by using the function A(T), 

^(-)^ ^5 («) 
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Figure 6: Two critical bubbles. Dotted line is the extremum solution for A = 0.9 
and solid line for A = 0.6. The dimensionless radial variable is r' = Mr, where M 
stands for the bosonic mass in the potential, M^{T) = j{T^-T^). 



which satisfies A(To) = 0, A(Tc) = 1. The nucleation action is written in terms of the 
function A(T) as follows: 

^(M^)) - ^^2jY^2 ■ (44) 

The function /(A) has been determined in Papers I and III by solving numerically the 
field equation (|12D for the potential (|39|). It is a smoothly behaving function with the 
special value /(I) = 1. In fig. 6, the bubble solution of the field equation is shown for 
two illustrative values of A. For A = 0.9 the solution somewhat resembles a thin- walled 
bubble, but for A = 0.6 the bubble core is far away from the true vacuum. 

From eq. (^) one can solve for tf, the age of the Universe when the bubbles had 
filled the space. Explicitly, the equation determining the amount of supercooling is, in 
analogy with demanding that the right side of eq. (0) equals unity, as follows: 

1 d^ 



A- S{Xf) -4\og 



, (45) 



_(T,/To)2-ldA 

where A is given by eq. (|3^) and Xf = A(Tj). The thin-wall limit, or small relative 
supercooling scenario, considered in Subsection follows as a special case if the 
solution A/ of the above equation is close to unity. For example, the expression (|3TD 
for the nucleation action follows from eq. (^) in the limit A/— >1. 
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Let us now inspect in more detail how the cosmological electroweak phase tran- 
sition is assumed to have proceeded. The sequence of events is illustrated in fig. 7. 
The numerical values of various quantities presented in the figure correspond to the 
following values of the parameters: 

= 100 GeV, 7 = 0.1309, 

(46) 

a = 0.0162, A = 0.0131. 

The values of 7, a and A are those obtained by Huet et al. [1993] by substituting 
the zero-temperature masses Mt = Mw, M/i = 40 GeV Mw/2 into the perturbatively 
evaluated effective potential. Here the unrealistically low mass for the Higgs-particle is 
more crucial than that for the top-quark, as will be discussed shortly. The parameter 
values (^6]), used here just for illustration, imply that the phase transition is very 
weakly first order. This is indicated by the small latent heat (= 0.086 T^'^) and the long 
correlation length (= IS.O/T^ in both phases at T^). 

Trust in the perturbative calculations leads to complications here. The Higgs mass 
used among others by Huet et al. [1993] is close to the upper limit of the allowed region if 
one requires that the baryon asymmetry of the Universe was created in the electroweak 
transition. A necessary condition for the electroweak baryon asymmetry generation is 
that the sphaleron transitions should have been frozen out in the low-temperature phase 
already at Tc. For that to have happened the Higgs scalar must not be heavier than 
about 40 GeV according to perturbative analysis [e.g., Huet et al. 1993]. On the other 
hand, the Higgs mass of such a low value seems to be ruled out by LEP experiments 
[Davier 1992]. As a solution to this dilemma some authors have assumed the existence 
of additional scalar fields. In multi-Higgs cases the parameter Mh is not the true zero 
temperature mass of a physical Higgs particle and is hence not constrained by the mass 
measurements [see Liu, McLerran and Turok 1992, and references therein]. 

In our phenomenological approach the parameter values quoted in eq. (^) do 
not pose any problems. Even within the minimal standard model these values do not 
imply a Higgs mass which were ruled out experimentally, as long as one does not try 
to extend the vahdity of the potential down to zero temperatures. 
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Figure 7: Schematic 1+1 dimensional figure of the electroweak phase transition. 
Grey area presents the high-temperature and white the low-temperature phase, 
and the separate pictures at right give two-dimensional snapshots of the transi- 
tion. The numbers are for the parameter values in eq. (|^). The idea for this 
figure is borrowed from Rummukainen [1990]. 



In fig. 7, three essentially different time scales can be seen. One is the Hubble 
time at the transition. Another is the duration of time the Universe stayed in the 
supercooled metastable state, Atgc = tf — tc- The phase transition time tf can be 



determined from eq. (45) by utilizing for the nucleation action the numerical function 
given in Paper III. At the nucleation, the value of the action turns out to be only 106 
instead of the naive value which according to eq. (|^) equals A (= 145), i.e., the last 
term in that equation does have some significance. The shortest time scale appearing 
in fig. 7 is the growth time of the shocks, defined as Atsh = -Rnuci/^- The average 
distance of nucleation centers -Rnuci is defined in eq. (^61). In numerical calculations the 
dimensionless prefactor b of the nucleation rate in eq. ([l^ was taken to be unity and 
the free-gas value v = l/ y/S was used for the shock velocity. 

The phase transition was completed at time te when the actual bubbles of new 
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phase, expanding as deflagrations behind the shocks, had met and coalesced fully. 
(Detonations, the other type of explosive processes, would require much more super- 
cooling than eq. ( ^5] ) indicates for the present parameters.) The growth time of the 
bubbles of the new phase cannot be solved, because the velocity of the deflagration 
front is not known. However, it seems probable that the difference of the velocities 
of the deflagration and shock fronts was clearly less than one order of magnitude 
[Ignatius, Kajantie, Kurki-Suonio and Laine 1993], which indicates that the growth 
time of deflagrations did not differ drastically from Atsh- This also means that the 



inaccuracies in the estimates of the time scales are not signiflcant (c/. Subsection [4.2|) . 

Mutual collisions of the shocks reheat the system to a certain temperature T^h. 
As will be demonstrated below, the reheating temperature Trh was in the electroweak 
phase transition less than the critical temperature Tc. Because of this and the fastness 
of the phase transition, Trh can be estimated by going to the extreme case where the 
whole space were converted instantaneously from the old to the new phase at time tf. 
DeGrand and Kajantie [1984] called this scenario "abrupt transition". The reheating 
temperature can be obtained from the equation 

e,(Trh) = e.iTj) . (47) 

For the parameter values given in eq. ( ^6] ) the reheating turns out to be only 13% on 
the scale where 100% would mean reheating back to the critical temperature, i The 
fact that the critical temperature is not reached during the reheating is crucial for 
the scenario of baryon asymmetry generation discussed by Liu, McLerran and Turok 
[1992]. If the reheating temperature Trh had been close to Tc, the velocity of deflagration 
fronts would have decreased substantially (see Subsection ^4.2|) . In that case the baryon 
number produced in the front would have had time to diffuse to the old phase, where 
it would have been washed out by the sphalerons. 

At the flnal stage of the phase transition the shrinking droplets of the old phase 
produced rarefaction waves; however, their effect was shadowed by the presence of the 
remnants of shocks. Recently Huet et al. [1993] have demonstrated that the deflagra- 



^In condensed matter physics a first-order phase transition with less than 100% reheating is called 
"hypercooled" instead of supercooled [see for example Leggett and Yip 1990]. 
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Figure 8: Relation between temperature and time in the electroweak (EW) and 
in the quark-hadron (QH) phase transition. Dotted curve denotes the reheating 
period during which temperature was far from being spatially uniform. 



tion fronts did not develop any instabilities while expanding. Thus the length scale of 
inhomogeneities is given by the usual expression (p6|) for the average distance of nu- 
cleation centers -Rnuci- However, it is hardly probable that the density inhomogeneities 
produced in the electroweak phase transition could have been of any importance. For 
instance, for the parameters ( ^6]) one can estimate the maximal relative pressure dif- 
ferences to be only Ap/p ^ 4 x 10~^. In comparison, for the storms in the atmosphere 
of Earth this quantity can be three orders of magnitude larger. 

In fig. 8 the behavior of temperature versus time both in the electroweak and in 
the quark-hadron transition is shown schematically in a log-log plot. From the figure 
we can see that the electroweak phase transition did not last long. Afterwards, the 
usual relation between temperature and time, given in eq. (^), became soon valid again. 
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4.2 The Quark— Hadron Phase Transition 



Investigations of the cosmological quark-hadron phase transition started over a decade 
ago [Ohve 1981; Suhonen 1982]. After these early studies it was soon reahzed that 
the onset of the supposedly first-order transition required supercooling [Hogan 1983]. 
Although much progress was made during the subsequent years in understand- 
ing various features of the transition [Witten 1984; DeGrand and Kajantie 1984; 
Applegate and Hogan 1985; Kajantie and Kurki-Suonio 1986; Fuller, Mathews and Al- 
cock 1988], several questions are still unanswered, partially due to insufficient knowl- 
edge of the properties of thermal quantum chromodynamics. (For a more complete 
reference list see Bonometto and Pantano [1993].) 

Lattice Monte Carlo simulations provide the best tool currently available for the 
study of the equation of state in QCD near the transition point. But even with this 
method one has not been able to solve the order of the transition for the cosmologically 
relevant case, i.e., when the chemical potential vanishes and there are two light (u,d) 
and one intermediate-mass (s) quark species. What is known is that in the case of pure 
glue, corresponding to infinitely heavy quarks, the transition is first order, and so it 
is with four light quarks. (For a review on lattice results see [Petersson 1992].) Also 
is known that there is a substantial jump in the energy density within a temperature 
interval of less than 10 MeV around the critical temperature. However, for working out 
the consequences of the phase transition in cosmology this information is not sufficient 
as long as one is not able to distinguish between a first-order, and a second-order or 
non-existing transition. If the transition is not first order, no supercooling can occur, 
even if the equation of state gave rise to a very rapid change in the energy density. 
This is due to the extremely slow expansion of the Universe. 

Guided by the recent lattice calculations, we use in numerical estimates the fol- 
lowing values for the physical quantities of the quark-hadron phase transition: 

Te = 150 MeV , L = 2T^ , a = 0.02 . (48) 

Here is the critical temperature, L the latent heat and a the interface tension. The 
true value of the critical temperature lies very probably somewhere between 100 and 
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250 MeV, and 150 MeV may be a good guess for its value [Petersson 1992]. For the 
other two quantities there is currently no lower limit, since they vanish if the phase 
transition is not first order. The values of L and a given in eq. ( ^8] ) are based on pure 
glue lattice simulations: the value of the latent heat is taken from Iwasaki et al. [1992], 
and the interface tension has been determined in computer studies by Grossmann and 
Laursen [1993], where the length of the lattice in the imaginary time direction was 2 
lattice points, and by Iwasaki et al. [1993], where it was 4 and 6 lattice points. 

The simplest analytical model for the QCD equation of state is that of the MIT 
bag model, which is often employed in the cosmological context. In this model the 
pressures of quark and hadron phases are given by 

p,{T) = a,T* - B , ph{T) = auT^ . (49) 

The value of the bag constant B is determined from the condition of equal pressure 
at Tci B = {ag — ah)T^. The bag equations (^) also follow as a special case from the 
more general equation of state presented in Subsection [4.1| . Approximating B{T) = 
(L/4)(l — T'^/T^), which in the limit of small supercooling is equivalent to eq. (|37D , 
and substituting this to the equation of state (^T]) with the constant —B{0) added in 
both sides, one recovers the bag equations: 

p,{T) = aT'-^, PK{T) = {a-^)T\ (50) 

The equation of state of the naive bag model, which follows from counting the 
particle species and utilizing eq. (|), gives an approximate upper limit for the latent 
heat. Somewhat above the critical temperature the strongly interacting relativistic 
particles are the gluons, and u- and d-quark. (Depending on its mass, also the s- 
quark could be counted.) For temperatures somewhat below Tc the only strongly 
interacting particles that we include within this naive approach are the pions, the 
lightest hadrons. The other low-massed and massless particles — the photon, electron, 
muon and neutrinos — are present in both phases. The effective number of relativistic 
degrees of freedom in the quark-gluon and in the hadron phase is thus given by g^^q = 
51.25 and 5'*ft = 17.25, respectively. 

In fig. 9, the energy density of the cosmic fluid is schematically plotted both for the 
naive bag model, and for a more realistic equation of state consistent with the lattice 
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Figure 9: Behavior of energy density in QCD. Thin curve is for the naive bag 
equation of state, and thick curve for a weaker first-order transition. Dotted hues 
denote the metastable branches. 

simulations. We clearly see how the naive bag model exaggerates the value of latent 
heat, Lbag = 14:.9T^ ^ L. The parameter values of the bag model can be corrected 
to reproduce a desired latent heat. However, the corrected bag model does not mimic 
well the realistic equation of state over the whole range. But in the vicinity of it 
reproduces the true equation of state with a first-order accuracy in the pressure and 
zeroth-order accuracy in the energy density, which is sufficient for determining the onset 
of nucleation. In this case the use of the bag model with improved parameter values 
is in practice equivalent to employing directly the thin-wall approximation discussed 
in Subsection |3.3| , except that changing the equation of state also affects the relation 
between time and temperature of the Universe. 

Validity of the thin-wall approximation is violated if the correlation length asso- 
ciated with the transition is not clearly smaller than the radius of the critical bubble. 
It is not easy to tell what the relevant correlation length is in QCD near the transition 
temperature. However, since the transition is just weakly first order, if first order at 
all, it is quite possible that the correlation length is quite large. 

In a case where the thin-wall approximation is inapplicable, one could employ in 
nucleation calculations the order parameter model presented in Subsection |4.1| . There 
is a one-to-one correspondence between the four parameters of the potential (|39D of 
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that model and the four physical parameters of the transition [Paper I; Kajantie 1992]: 



[To, 7, a, A] < — > [Tc, a, L, y , (51) 

where Ic is correlation length. (In the order parameter model with a quartic potential 
the correlation lengths in both phases are equal to Ic at the critical temperature.) Once 
the values of the physical parameters are known, also the parameters of the potential 
are completely fixed. 

It is not clear how one should interpret the order parameter field in QCD because 
the theory does not have any classical potential driving the transition. If the latent 
heat is not small, the order parameter represents several degrees of freedom. Then it 
would seem more natural to identify the 0-field with a thermodynamical quantity like 
the energy density, in the same manner as was done by Csernai and Kapusta [1992a]. 

From now on we will assume that the nucleation in the cosmological quark-hadron 
phase transition took place under conditions which were close to the thin- wall limit. 
For the values of the physical parameters presented in eq. (^8|) , eq. (|36| ) gives for the 
critical bubble a very large value of the radius, Rcr{Tf) = 38/Tc = 51 fm. This shows 
that the thin-wall approximation would remain valid even if the correlation length were 
large. 

The main events of the quark-hadron phase transition are shown in fig. 10 in 
the same way as was done in fig. 7 for the electroweak case. Inspecting first the 
early stages of the phase transition, we note that in the quark-hadron phase transition 
the supercooling was smaller than in the electroweak transition. Secondly, we may 
compare the growth time of shocks with the duration of supercooling utilizing for 
example eqs. (0) and (|35|) : 

— ^ = ■ 52 

Atsc A ^ ' 

It is interesting to note that this ratio is completely determined by the nucleation action, 
or approximately by the age of the Universe. In other words, at any temperature there 
is the definite relation (^2]) between the degree of supercooling and the nucleation 
distance (assuming that the value of the shock velocity is a constant). This result is in 
principle valid only in the thin- wall limit. However, from fig. 7 one can infer that the 



41 



prediction holds rather well also for the electroweak transition in the example case we 
considered. 

The values derived for the two time scales Atsc and Atsh may be somewhat er- 
roneous, since in three dimensions the shocks are weak, especially if the deflagration 
front is very slow [Kurki-Suonio 1985]. It seems that the deflagration front velocity 
was indeed quite small, probably at least one order of magnitude smaller than the 
velocity of the shock front [Kajantie 1992; Ignatius, Kajantie, Kurki-Suonio and Laine 
1993]. The weakness of the shocks could make the estimates of the nucleation process 
inaccurate, because new bubbles could possibly nucleate to a region already touched 
by a shock. This inaccuracy can be at least partially cured by using in the calculations 
an effective velocity v, which is smaller than the true velocity of the deflagration front. 

The main difference between the electroweak and the quark-hadron transition 
is that only in the latter transition the Universe was reheated back to the critical 
temperature (see flg. 8). This is due to the much larger value of the latent heat in 
the quark-hadron transition. After the reheating, the phase transition proceeded very 
slowly, and almost in thermal equilibrium. The expansion of the Universe did not 
cause any cooling; instead, the denser quark-gluon matter was transformed to the 
more dilute hadron matter. As is discussed in Paper I, the duration of this period can 
be approximately determined from the relation 

fr^" • ^^^^ 

This approximation holds if the resulting value is clearly smaller than unity. In the 
case of this period of slow burning the expansion rate of the Universe determines 
the typical velocities of deflagration fronts, too. The velocities are roughly given by 
-Rnuci/-Rhor divided by L/Asq. The numerical value that this gives for the velocity is of 
the order of 10""^. 

At the flnal stages of the phase transition the decaying quark droplets produced 
rarefaction waves [Kajantie and Kurki-Suonio 1986]. This led to the creation of density 
inhomogeneities, which in principle could have signiflcantly affected the nucleosynthesis 
and could be observed in the present-day Universe. However, if the parameter values in 
eq. (|48|) are roughly correct, the distance scale of density inhomogeneities was too short 
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Figure 10: Schematic 1+1 dimensional figure of the phase transition from the 
quark-gluon (grey) to the hadron phase (white). This figure was presented orig- 
inally by Rummukainen [1990]; the current version is a modified one. For clarity 
only the effect of the dying quark droplets is shown in the world lines of the cos- 
mic fluid (thin dotted lines), the effect of an expanding hadron bubble has been 
illustrated earlier in fig. 3. The numbers correspond to the values of physical 
quantities in eq. (48). 
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for this to happen. Only a distance scale -Rnuci of at least one meter at the quark-hadron 
transition temperature could have had later an effect on the abundance of light elements 
[Applegate, Hogan and Scherrer 1987; Kurki-Suonio et al. 1990; Mathews et al. 1990]. 
Redshifted to the present Universe, the length -Rnuci ~ 4 cm is less than the distance 
from Earth to Sun. In the cosmic scale this is a very short distance. 

If the interface tension were the same as in eq. ( ^Sf ) but the latent heat much 
smaller, the distance scale -Rnucb given in eq. ( ^5|) and proportional to a^^'^/L, would 
correspondingly be much larger. This possibility cannot presently be ruled out, because 
the values of the physical quantities a and L are not known. However, it is tempting to 
think that the values of these quantities would not be arbitrary if the phase transition 
were very weakly first order; that instead in this limit they would show some sort of 
universal behavior compared with other physical transitions. By weakly first order 
transition we mean a transition in which the values of the thermodynamical quantities 
L, a and 1/1^ are in dimensionless units small when compared with unity. 

A phase transition which can be treated analytically is the one between ordered 
and disordered phases in the two-dimensional g-state Potts model [see e.g. Wu 1982]. 
This model is a generalization of the Ising model to q spins, and the transition is 
for q > 4 first order, the stronger the larger q is. In the case of the transition 
of the two-dimensional Potts model there are no additional parameters besides q. 
The values of the latent heat [e.g., Wu 1982] and more recently, the interface ten- 
sion [Borgs and Janke 1992], have been calculated analytically. In the limit where the 
phase transition becomes weaker and weaker, that is, when q approaches 4 from above, 
the interface tension vanishes more rapidly than the latent heat. (A similar behavior 



can be seen in the quartic order parameter model presented in Subsection [4.1| , if one be- 
lieves in the naive way of weakening the phase transition: by decreasing the value of a 
and keeping the other parameters appearing in the potential (^9]) constant one observes 
that cr^/^/L oc a ' .) The scaling argument coming from the two-dimensional Potts 
model seems to hint that if a cosmological first-order phase transition is made weaker, 
the distance scale between nucleation sites gets smaller — a quite natural behavior. 

It seems that in the cosmological quark-hadron phase transition the expand- 
ing deflagration bubbles were on the borderline between stability and instability 



44 



[Huet et al. 1993], if one assumes that latent heat was carried away from the front 
by hydrodynamic flow and not by neutrinos. The opposite assumption has also been 
made [Preese and Adams 1990; Adams, Preese and Langer 1993], and it easily leads to 
instabilities. However, according to Applegate and Hogan [1985] the relative impor- 
tance of neutrinos as heat carriers vanishes in the limit of small supercooling. 

Finally, let us mention two exotic topics related to the cosmological quark-hadron 
phase transition. Firstly, Mardor and Svetitsky [1991] made the observation that in 
the MIT bag model small hadron bubbles exist already above Tc. But this does not 
imply that the same would be true for the real QCD. Indeed, employing the quartic 
order parameter model the normal behavior is recovered (Paper III). Secondly, in some 
circumstances it might be possible for the quark droplets to survive over the transition 
[Wittcn 1984]. If stable, these lumps of strange quark matter could then in principle 
be observed in the present-day Universe. 
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5 Conclusions and Future Prospects 



The goal of this Thesis has been to achieve an insight to the physical processes that 
occurred in the Universe during different first-order phase transitions. 

We have pointed out that especially the onset of a cosmological phase transition 
shows a universal behavior allowing for a general approach. We have presented in this 
work some general methods which can be applied to an analysis of the sequence of 
events taking place in cosmological first-order phase transitions. 

In the case of the cosmological electroweak or quark-hadron phase transition, the 

quantitative description has not yet reached a fully satisfactory level. The reason for 
this is the lacking knowledge of the correct input values of the physical parameters. 
The case with the phase transition suggested in Paper II is different in the sense that 
there the main uncertainty comes from the basic assumptions of the scenario itself. 

In the future one will hopefully learn the correct values of the physical quantities 
related to the electroweak and quark-hadron phase transitions. With their improving 
accuracy lattice simulations should provide this information for QCD with physical 
quark masses, as well as for the electroweak theory. 

Combining lattice simulations with analytical calculations in effective three- 
dimensional theories, it seems to be possible to derive for the electroweak the- 
ory a potential incorporating even non-perturbative effects [Shaposhnikov 1993; 
Farakos et al. 1993]. Such a potential could then be employed for studying the electro- 
weak phase transition. 

In order to understand in detail the dynamics of cosmological phase transitions 
one must have a good knowledge of the bubble growth at microscopic level. At present 
we are investigating the growth of bubbles using a model in which there is a friction- like 
coupling between the order parameter field and a cosmic fluid field. By employing this 
model one is able to numerically simulate for instance the collisions between bubbles. 
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Moreover, the velocity of a deflagration wall can be determined exactly as a function of 
the friction coefficient. So far, we have applied the model in 1+1 dimensions [Ignatius, 
Kajantie, Kurki-Suonio and Laine 1993]. In the future, the computations should be 
extended to include spherically symmetric three-dimensional bubbles. Furthermore, 
one could try to obtain a good estimate for the value of the friction coupling that 
determines the velocity of the bubble wall. 

By combining all these developments, one is led to the conclusion that it is pos- 
sible to achieve an accurate quantitative description of different cosmological phase 
transitions in the not too distant future. It would then also be possible to calculate 
the value of the baryon asymmetry of the Universe from the first principles. 
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